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Abstract. We present BeppoSAX X-ray spectra of the rem- 
nant of the supernova of AD 1006, which cover a broad spec- 
tral range of 0.1-10 keV. In our analysis we concentrate on the 
thermal emission from the central region of the remnant. For 
this purpose we fitted the spectra of the whole remnant using 
a new version of the spectral program SPEX, which takes into 
account the interdependence of spectra of different regions as 
the result of the point spread function of the instruments. This is 
in particular important for separating the synchrotron emission, 
coming from the bright rims of the remnant, from the thermal 
emission from the central regions. 

The thermal emission appears best fitted with a combina- 
tion of two spectral components with electron temperatures of 
^0.7 keV and >, 3 keV. The plasma in SN 1006 appears to be 
very far out of ionization equilibrium (n c t ~ 2 10 15 m~ 3 s), 
which means that the L-shell emission of magnesium, sil- 
icon and sulphur contribute significantly to the flux below 
0.5 keV. We also note the presence of Fe K line emission at 
6.3 ± 0.2 keV, in agreement with an origin from inner shell ion- 
izations and excitations. We confirm that the abundances in SN 
1006 are clearly non-solar. In particular silicon is very abun- 
dant. We did not find substantial variations in plasma proper- 
ties over the face of the remnant, although the hotter component 
seems more dominant in the Northern half of the remnant. 

Key words: ISM: individual objects: SN 1006 - ISM: super- 
nova remnants - X-rays: ISM 



1. Introduction 

The historical supernova SN 1 006 w as probab ly a T ype la su- 
pernova (Clark & Stephenson 1977 , Schaefer 1996 ). Its rem- 
nant has now a diameter of 30' and in recent years three inter- 
esting findings were reported: 

- The ASCA X-ray spectrum indicates that the emission 
above ~ 1 ke V is dominated by synchrotron radiation from 
electrons with energies up to ^100 TeV, which are accel- 
erated at the shock front (Koyama et al. 



the presence of extremely relativistic electrons (Tanimori 
et al. 



1998). 



1995, Reynolds 



1998). The detection of TeV gamma ray emission confirms 
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- Modeling far ultra-violet spectra of a Northwestern 
filament has revealed that the post shock plasma has 
not reached ion-electron temperature equilibration. The 
electron temperature was f ound to be only 10% of the ion 
temperature (Laming et al. 



1996) 



- Ultraviolet absorption feature s in the spectrum of the 
Schweizer-Middleditch ( 1980 ) star, a blue subdwarf lying 
behind the remnant, have revealed the presence of blue 
and red-shifted unshocked silicon and iron (Hamilton et al. 



1997 ). Remarkably, no evidence was found for the pres- 
ence of ~ 0.5 Mq of iron, which is likely to have been syn- 
thesized during the carbo n defl agration of the white dwarf 
progenitor (Nomoto et al. 1984 ). Carbon deflagration is the 
currently favored model for Type la supernovae. 

The presence of a synchrotron component contaminates the 
thermal X-ray emission and the non-equilibration of ion and 
electron temperatures obscures the underlying shock hydrody- 
namics, which might otherwise have been inferred from the 
temperature structure. On the other hand, in X-rays other ion 
species can be studied than in the UV and this may help veri- 
fying the models for nucleosynthesis of Type la supernovae. 

In this paper we present BeppoSAX spectra of the rem- 
nant of SN 1006. We will show that the thermal emission from 
the central region of the remnant is best described by a two 
component non-equilibrium ionization model with a very low 
ionization parameter. 

2. The data and method 



The Italian-Dutch satellite BeppoSAX (Boella et al. 1997a) ob 



served SN 1006 in April 1997. The observation consists of 
three different pointings, which together cover the whole rem- 
nant by the imaging instruments LECS (Parmar et al. 1997) 



and MECS (Boella et al. |1997b[ ). The total observation time is 
91 ks for the MECS and 37 ks for LECS, which has a lower 
observation efficiency owing to UV leakage. 

Both the LECS (~ 0.1 - 10 keV) and the MECS (~ 
2 — 10 keV) are gas scintillation counters with a spectral reso- 
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Fig. 1. The right images shows the sky regions for the model spectra. The image on the left shows a ROSAT HRI mosaic of SN 
1006 on the same scale. 



lution of ~ 8% at 6 keV, a spatial resolution at 6 keV of ~ 2' 
FWHM and a half power width]] of ~ 2.5'. The spatial res- 
olution is strongly energy dependent: with increasing energy 
the core of the point spread function decreases, whereas the 
wings of the point spread function (determined mostly by the 
mirror properties) become more dominant (Parmar et al. 1997, 
Boella et al. |l997b[ cf. Vink et al. |1999| ). The extended wings 
of the point spread function make that in the case of SN 1006 
the spectra of the low brightness central region, which emits 
mostly thermal emission (Koyama et al. 1995, Willingale et al. 



1996, the latter based on ROSAT PSPC data), is contaminated 



by the radiation from the rims of the remnant, which is pre- 
dominantly synchrotron radiation. 

We circumvented this problem by using a new version 
(v2.0) of the X-ray spectral fitting program SPEX (Kaastra et 
al. 1996, cf. Kaastra et al. 1999), with which one can fit several 



spectra from the same object simultaneously. In this version the 
concept of the spectral redistribution matrix is extended to in- 
clude also a spatial redistribution part. The spectral model is 
calculated on a spatial/spectral input grid and a model of the 
instrument properties describes what fraction of the incoming 
photons with a certain energy, coming from a certain region of 
the sky will end up in a given energy bin (channel) and spa- 
tial bin (spectral extraction region). This means in practice that 
two sets of spatial regions have to be provided (in the program 
these are specified using images compliant with the FITS for- 
mat). One set specifies the sectors for which the model is calcu- 
lated, the other set specifies the spectral extraction regions. For 



The half power width gives the width of the area in which 50% of 
the photons of a point source are expected to end up; this measure is 
sensitive to the broad wings of the point spread function, whereas the 
FWHM is more determined by the core of the point spread function. 



stable calculations, the extraction regions should roughly cor- 
respond to the spatial model sectors, but a detailed correspon- 
dence is not necessary. For example, some regions of SN 1006 
were only partially covered by some observations, but this is 
accounted for in the spectral/spatial redistribution matrix. The 
X-ray emission models used for this paper are the same as in 
the standard SPEX program. 

We divided the remnant into six regions. Apart from a cen- 
tral region and two X-ray bright rims, we divided these regions 
further into a Northern and Southern half in order to see if 
the thermal emission from the brighter Southern half is dif- 
ferent from the emission from the Northern halve. With North 
and South we mean here, and throughout the rest of the text, 
North and South with respect to the principle axis of the rem- 
nant, which is tilted with respect to the North-South direction 
in equatorial coordinates. The spatial model grid for which the 
models were calculated are shown in Fig. [j] We chose spec- 
tral extraction regions roughly corresponding to the model re- 
gions, but the extraction regions corresponding to the bright 
rims were larger, in order to reduce the contamination of the 
central regions by emission coming from the bright rims. We 
used archival ROSAT HRI images (Fig. [j], see also Winkler & 
Long 1997) to define the spectral extraction regions and the 
model sectors. For the spectral extraction regions we took also 
into account the extend of the synchrotron rims in the Bep- 
poSAX images, as they seem broader than on the HRI images 
due to the point spread function. 

The extraction of the spectra and the generation of the 
spectral-spatial response matrix was done with a program 
specifically made to extract BeppoSAX spectra. It incorporates 
the current knowledge of the detectors, such as the point spread 
function, vignetting and absorption by the strong-backs. We 
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Fig. 2. The combined LECS/MECS spectrum of the Southern 
central region, with a simple model fit (i.e. no hard thermal 
or non-thermal component), which does not adequately fit the 
spectrum above ~ 5 keV. 

used the response matrices of September 1997 and the standard 
background data (November 1998). However, there are still 
some uncertainties in the detector calibration, especially for 
off axis positions. These uncertainties are larger for the LECS, 
which has a more complicated design than the MECS. For this 
reason we left the normalization of the LECS spectra and some 
off-axis MECS spectra free with respect to the on-axis MECS 
spectra. The relative LECS normalizations turned out to be be- 
tween 0.6 to 0.9, co nsistent with other BeppoSAX results (e.g. 
Favata et al. 1997). The spectra were binned to a bin size of 
roughly l/3rd the spectral resolution and some further rebin- 
ning was done for channels with low count rates. In order to 
circumvent statistical problems with bins with very few counts 
we used a method proposed by Wheaton et al. (1995). This 
means that after obtaining a good fit, we used the best fit model 
to calculate the expected error per bin, instead of the observed 
counts. Using this method with two or three extra iterations 
gives in general a stable and in principle more reliable \ 2 value. 

3. Spectral fitting 

Fitting six spectral regions simultaneously has the obvious dis- 
advantage that the spectral model can become very complex. 
So we took care to constrain the spectral model as far as pos- 
sible without loosing too much of its heuristic qualities. For 
each sector we chose to have three to four spectral components. 
A typical configuration consisted of the following components 
for each sky region: a power law component, one or two non- 
equilibrium ionization (NEI) thermal components and an ab- 



sorption component (Morrison & McCammon 1983). We only 



looked for spectral differences in the thermal emission between 
the Northern and Southern halves of the remnant. We assumed 
uniform abundances for SN 1006. 

Based on their deprojection of th e ROSAT PSPC image 
of SN 1006, Willingale et al. ( |l996| ) reported that the syn- 
chrotron emission does not seem to originate from all around 



Table 1. The best fitting parameters for our model with two 
NEI thermal components. The power law components for the 
emission from the rims are listed separately in Table || The 
abundances are relative to solar abundances (Anders & Grevese 



1989). Errors are 90% confidence limits and lower/upper limits 



are 99% confidence limits. 



North 



South 



n e nnV/d 2 I (10 62 m- 3 kpc -2 ) 

T c I (keV) 

n c t I(10 15 m~ 3 s) 

n c n n V/d 2 II (10 62 m~ 3 kpc^ 2 ) 

T e II (keV) 

n c t II (10 15 m" 3 s) 

C 

O 

Nc 

Mg 

Si 

s 

Fc 

JV H (10 20 cm" 2 ) 



2.3 ± 1.0 
0.71 ±0.15 

3.3 « 4.6) 

0.1 ±0.04 

3.4 (> 3) 
0.7 « 2) 



3.9 ± 1.0 
0.78 ±0.09 
2.3 « 3.3) 

0.14 ±0.03 
8.0 (> 4) 
0.7 «3) 



0.8 ±0.4 
0.34 ±0.07 
0.28 ± 0.09 
0.8t°;» 
9.1 ± 1.5 
5 ± 2 
4±2 

8.8 ±0.5 

1678/1285 



the remnant, but only from incomplete shells at the Northwest 
and Southeast of the remnant. Consequently, in our simplest 
model we assume that no synchrotron emission is originating 
from the central region of the remnant. This model produces 
a reasonable fit to the LECS spectra, but it does not fit the 
MECS spectra of the central regions above ~4 keV, where an 
excess in the observed spectra with respect to the model ex- 
ists (Fig. fy. The temperature of the thermal components was 
T c ^ 1.5 keV. So clearly an additional emission component is 
needed in order to fit the emission from the central regions. 
This can be either a thermal component with an higher tem- 
perature, but it could also mean that the non-thermal emission 
observed to come from the rims has in reality cylindrical sym- 
metry, in which case the apparent structure of the remnant in 
X-rays may be due to an extreme case of limb brightening. 
The latter possibility would be in disagreement with the above 
mentioned ROSAT PSPC findings, but it is conceivable that the 
synchrotron emission is coming from suc h a thi n layer that the 
deprojection scheme of Willingale et al. (1996 1 may not have 
worked adequately. 

We investigated both possibilities. In the case of an addi- 
tional power law component we fixed the power law index of 
the central regions to a value of 2.8, similar to the the values 
found for the rims (see Table ||). We found, however, that an 
additional thermal component offers a better explanation. Only 
a hot thermal component fits adequately the Fe K emission seen 
in Fig. [3] (most clearly in the right bottom panel). This emission 
implies that there is some plasma present with a temperature in 
excess of ~ 2 keV. An additional power law component has, 
moreover, the disadvantage that the result is at a closer look 



4 



Jacco Vink et al.: The BeppoSAX X-ray spectrum of the remnant of SN 1006 



Table 2. The best fit parameters for the non-thermal emission 
of SN 1006. 





norm 


index 




(s _1 m _2 keV -1 @ 1 keV) 




Northeast 


181 ± 10 


2.81 ±0.05 


Southeast 


83 ± 10 


2.76 ±0.08 


Northwest 


86 ± 13 


2.90 ±0.16 


Southwest 


76 ±9 


2.83 ±0.12 



inconsistent: it results in very high abundances for the thermal 
component (e.g. 148 times solar for silicon and 356 times solar 
for iron), whereas the emission measure indicated a total mass 

Of' 



5 M Q for a distance of 1.8 kpc (Laming et al. 1996). This 
implies that most of the shock heated plasma is swept up ma- 
terial, if SN 1006 was a Type la with a mass of ~ 1.4M Q , but 
then it is not clear why the abundances are so high. High abun- 
dances are also hard to reconcile with the evolved dynamical 
status of the remnant, which implies that the remn ant co nsists 
most ly of s wept up interstellar matter (Moffet et al. 1993 , Long 
et al. 1988). However, X-ray synchrotron models for supernova 



remnants (Reynolds, 1998 ) show that the synchrotron compo- 
nent is likely to have a more curved spectrum, which predicts 
less flux at lower energies. It may well be that such more ad- 
vanced models do not have the problems mentioned above, al- 
though an additional thermal component will still be needed for 
producing the Fe K emission. From a statistical point of view 
the two models fit the data equally well, with slightly better fits 
for the thermal model (the model with an additional power law 
has X 2 /v = 1689/1287). 

The best fit model with an additional thermal component 
(Table | and |) is shown in Fig. |[ The total emission measure 
of n e nuV/d 2 = (7.4 ± 1.4) 10 62 m~ 3 kpc~ 2 implies a mass 
of M = 8.3 ± 0.8/ 3/ 4 2 ^ 8 2 Mq, where f 0A i s the volume fill- 



ing factor divided by 0.4 (cf. Willingale et al. 1996) and di.g is 
the distance in units of 1.8 kpc (Laming et al. 1996| ). The vol- 
ume estimate assumes a spherical remnant with a diameter of 
30'. For this model the implied silicon mass is 0.05 M Q and the 
iron mass is 0.06 M Q , which is much lower than the 0.5 M Q 
of iron expected in remnants of Type la supernovae. The iron 
abundance reported here is, however, higher than reported by 



Koyama et al. (1995). We do not find significant variations in 
temperature or ionization between the Northern and Southern 
half of the remnant, but, as Fig. [3| shows, the hottest compo- 
nent seems more dominant in the Northern region. An inter- 
esting feature in the spectra of the central regions is the iron 
K-shell emission around 6.4 keV. This is especially apparent in 
the lower right panel of Fig. 0. The centroid of this emission 
(6.3 ± 0.2 keV) is consistent with emission from iron in low 
ionization stages (caused by inner shell ionizations and excita- 
tions). In the next section we will discuss the very low value of 
the ionization parameter (n e t ). 

Our best fit absorption colum n is a factor two lower than 
the value found by Koyama et al. ( 1995 ), but is higher than the 



value reported by Willingale et al. ( |1996[ ). The higher ASCA 
value is not very surprising, because ASCA was less reliably 
calibrated below ~ 1 keV and ASCA is not sensitive below 
0.4 keV. The fact that the ROSAT PSPC and BeppoSAX LECS 
show that there is emission at 0.2 keV is, however, a clear indi- 
cation that the absorption column is not as high as indicated by 
the ASCA spectrum. We think that our absorption estimate is 
consistent with the ROSAT PSPC spectra if one takes into ac- 
count that we use non-equilibrium ionization models and in this 
case such a model produces more line radiation below 1 keV 
than an equilibrium model, which implies a higher absorption 
to explain the flux levels. 

The power law index of the non-thermal spectra are some- 
what lower than indicated by the ASCA spectra and higher than 
found for the ROSAT PSPC data (see Table g). We checked 
for gradual steepening of the s pectru m, which is expected on 
theoretical grounds (Reynolds 1998 ). Indeed, we found some 
evidence that on average the power law index changes from 
2.65 ± 0.21 below 2 keV to 2.81 ± 0.05 above 2 keV. How- 
ever, the scatter in the four points is rather large and we can, 
furthermore, not exclude that the difference in slope is due to 
calibration uncertainties, as there are still problems with the in- 
tercalibration of the LECS and MECS instruments. 

4. Interpretation 

Our fits indicate that the spectra of the central region are best 
fitted with two NEI components with a very low value for the 
ionization parameter. The very low ionization is not surprising 
as previous studies indicate that the density of the ISM sur- 
rounding the remnant is very low: estimates of the pre-shock 
hydro gen number density vary from 0.04 cm -3 (Laming et al. 
1996 ) to 1 cm~ 3 (Winkler & Long 1997 ). Our mass estimate 
implies a pre-shock density of 0.1 cm~ 3 . For the Sedov model 
the average ionization parameter (n c t ) is approximately not, 
where no is the pre-shock hydrogen density and t the age of 
the remnant. So for the shocked interstellar gas we estimate n e t 
~ 3 10 15 m~ 3 s. For the ionization parameter of the shocked 
ejecta we can find an upper limit by assuming that 1 .4 M 9 of 
ejecta is completely ionized and that it has a volume filling fac- 
tor of 0.4. This gives an electron density of 0.04 cm~ 3 and 
implies n c t £ 1.2 10 15 m~ 3 s. The n c t value found by fitting 
the X-ray spectra agrees with our current understanding of SN 
1006. 

The low ionization has, however, some intriguing aspects. 
For example, our models indicate that most of the silicon is in 
the form of Si V to Si XII. So with the silicon lines around 
1.8 keV we are only observing the top of the iceberg. However, 
the silicon L-shell emission contributes significantly to the flux 
between 0. 1 ke V and 0.4 ke V. Also L-shell emission from mag- 
nesium, sulphur and argon, and K-shell emission from car- 
bon contribute to the flux below 0.5 keV. Unfortunately, L- 
shell emission is notoriously complicated; only the iron L-shell 
emission has been investigated in substantial detail. The LECS 
spectral resolution below 0.5 keV is not good enough to see 
possible discrepancies in the atomic data, but the resulting un- 
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certainties in the flux below 0.5 keV may have affected the 
spectral modeling. This demonstrates the significance of the 
broad spectral range covered by the two BeppoSAX instru- 
ments. The observation of emission from iron around 6.4 keV 
illustrates that even ions in low ionization stages emit some 
X-ray line emission, because of inner shell ionizations and ex- 
citations (in this respect SN1006 seems to be similar to RCW 



86, see Vink et al. 1997). It would be interesting to search for 



similar emission features from Ar and Ca with future X-ray 
detectors. 

The apparent two temperature structure is not surprising. 
SN 1006 is a large remnant and it would be more surprising if 
the whole remnant could be characterized by a single temper- 
ature. Furthermore, there are several possible explanations for 
the two temperature structure. Hydrod ynamical shock models, 
such as the Sedov model (Sedov 1959 ) predict temperatu re gra - 
dients. Models including a reverse shock (e.g. Chevalier 1982 ) 
even predict a two temperature regime. The lowest tempera- 
ture is associated with the ejecta heated by the reverse shock, 
whereas the highest temperatures are attained by the shocked 
interstellar medium. However, hydrodynamical models predict 
the ion temperature, whereas the shape of the thermal X-ray 
continuum is determined by the electron temperature, which, 
in absence of electron-ion equilibration, can in principle be as 
much as a factor 1835 (the electron/proton mass ratio) lower 
than the ion temperature. It can be even more if the plasma con- 
sist of highly enriched supernova ejecta. In these extreme cases, 
however, we would not observe thermal X-ray emission at all. 
The equilibration process is still poorly understood. In the case 
of SN 1006 there is strong evidence for poor equilibration of 



the electron and ion temperature (Laming et al. 1996). The ev 



idence comes from UV spectroscopy of the Northern shock 
front, but the equilibration may very well vary over the rem- 
nant. Another possible cause for temperature variations over 
the remnant are density differences. There is some evidence 
that SN 1006 has a very strong density gra dient f rom the front 
to the back of the remnant (Hamilton et al. 1997). Yet another 



interesting possibility is that the shock heating process does not 
produce a Maxwellian electron distribution, since the low den- 
sity in SN 1006 may preve nt a ra pid thermaliz ation o f the elec- 
tron distribution (Laming 1998, see also Itoh 1984, Hamilton 
& Sarazin 1984). The hot component could then be associated 
with bremsstrahlung from the tail of the non-thermal electron 
distribution. 

For practical reasons our fitting model has been largely a 
phenomenological model. It is important to realize this when 
interpreting the best fit parameters. For instance, our use of sin- 
gle/double NEI components does not account for the fact that 
there are temperature and ionization gradients. Or to be more 
precise: we suppose that these gradients can be approximated 
by two thermal plasma components. On the other hand mod- 
els incorporating temperature and ionization gradients, such 
as spectral models based on the Sedov model (e.g. Kaastra 



& Jansen 1993), suppose that the gradients are well defined, 



interstellar medium or because of the non-equilibration of ion 
and electron temperature, as discussed above. In our opinion a 
more severe problem with our modeling is that we use uniform 
abundances for the whole remnant, which is probably not cor- 
rect, as there is very likely freshly shocked ejecta present con- 
sisting entirely of metals, on the one hand, and shocked inter- 
stellar gas with solar abundances, on the other hand. This might 
for example affect our total mass estimate, which may be lower 
because a pure metal plasma is an efficient bremsstrahlung 
emitter. However, for more realistic models we need either data 
from more sensitive X-ray instruments, or we need a clear un- 
derstanding of the thermal and abundance structure of the rem- 
nant. For instance, with future X-ray missions we may be able 
to identify spatially or spectrally an ejecta emission component 
and observe the (projected) temperature structure of the rem- 
nant directly. 

The overall structure of SN 1006, as emerging from other 
observations, appears complicated and is not very well under- 
stood. Some observations seem even contradictory. For exam- 
ple, UV absorption measurements indicate that a lot of the 
ejecta are still in free expansion (Hamilton et al. 1997), but 



the kinematics of the remnant, on the other hand, show that SN 
1006 is dynamically evolved implying that most of the ejecta 



has been shocked by now (Moffet et al. 1993, Long et al. 1988). 



Finally, it is peculiar that a remnant that has such a high de- 
gree of symmetry with respect to the Southeast/Northwest axis 



(Roger et al. 1988 ) seems to have a strong front-back asymme- 



try (Hamilton et al. 1997). A strong deviation from c ylindri cal 



symmetry has also been reported by Willingale et al. ( |1996| ). 



5. Conclusion 

Our analysis indicates that a one temperature model does not 
fit the thermal spectrum of SN 1006 adequately. An additional 
thermal component with a temperature in excess of 3 keV is 
needed. The alternative, an additional power law component, 
cannot be excluded, but in this case the abundances of the low 
temperature thermal component become exceptionally high, 
which is inconsistent with the idea that SN 1006 must have 
swept up a considerable amount of interstellar material. 

The values for the ionization parameters for the thermal 
components appear to be very low (n c t ^ 3 10 15 m~ 3 s), but 
are consistent with the electron densities that we derive from 
the emission measure. Note, however, that in reality a range of 
ionization values will be present. For instance, the Sedov model 
(see Kaastra & Jansen 1993) predicts a superposition of plasma 



whereas in reality the gradients may very well be more com- 
plicated, because of the pre-supernova density structure of the 



ionization values ranging from roughly 0.2 in the center of the 
remnant up to 2 times the average value at a radius of 0.9 times 
the shock radius. Inhomogeneities in the pre-shock interstel- 
lar medium may extend this range upward and downward. For 
such a low ionization parameter the L-shell emission of Mg, 
Si and S contributes substantially to the flux below 0.5 keV. 
This explains that our best fit hydrogen absorption column (TVh 
= (8.8 ± 0.5) 10 20 cm ~ 2 ) turns out to be higher than found by 
Willingale et al. ( |l996j ). 
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Noteworthy is an iron K shell emission feature in the spec- 
trum. The centroid of the line emission is 6.3 ± 0.2 keV, which 
indicates that the emission is the result of inner shell ionizations 
and excitations. The fitted abundances clearly indicate non- 
solar abundances with, in particular, a very high silicon abun- 
dance. This is in agreement with the analysis of ASCA data 



(Koyama et al. 1995 1. We did not find significant differences in 



the temperatures and ionization parameters for the spectra of 
the Northwestern and Southeastern halves of the remnant, but 
the hottest component seems more dominant in the Northern 
region. 

Observations in the near future with Chandra (AXAF), 
XMM and Astro E will be able to reveal new information on the 
remnant which is needed in order to obtain a consistent model 
for this remnant. In particular, a spatial or a spectral separation 
of emission originating from shocked ejecta and shocked in- 
terstellar gas will be very useful for understanding more about 
the structure and abundances of this interesting but rather enig- 
matic supernova remnant. 
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Fig. 3. The spectra of various regions. The spectra from different pointings are here summed together for displaying reasons. The 
six upper panels show the spectra as observed by the LECS instrument and the six bottom panels are the spectra as observed by 
the MECS. The solid line shows the best fit model convolved with the instrument resnonse (see Table ET). The dotted line in the 
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